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AXIAL AND CIRCUMFERENTIAL VARIATIONS OF HOT-GAS-SIDE HEAT- 
TRANSFER RATES IN A HYDROGEN-OXYGEN ROCKET 
by Ralph L Schacht and Richard J. Quentmeyer 
Lewis Research Center 

SUMMARY 

An experimental investigation was conducted at the NASA Lewis Research Center to 
determine the axial and circumferential variations of heat -transfer coefficients in two 
rocket thrust chambers. Heat -transfer rates were determined from transient tempera- 
ture measurements at 20 locations in one thrust chamber and at 18 locations in another 
thrust chamber. Heat-flux meters were positioned at six circumferential locations in 
both the chamber and throat stations. The thrust chambers were operated over a range 

of mixture ratios from 2. 57 to 5.67 (28 to 15 percent hydrogen) at a nominal chamber 

o 

pressure of 2.068 MN/m (300 psia). Three injectors were used. 
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Data are correlated on the basis of two equations of the form St* Pr* ' = 

CRet - ®'^ and St* Pr*®' ^ = C Re* ~®'^, where St*, Pr* , and Re* are reference 
Stanton, Prandtl, and Reynolds numbers, respectively, and C is a constant. The length 
dimension used to determine the Reynolds number was either the station diameter d or 
the length from the injector face to the station axial location x. The maximum circum- 
ferential variation of C in the chamber and throat are 38 and 27 percent, respectively, 
when Re^ was used in the correlating equation. The mean values for C for the cham- 
ber and throat circumferential surveys were 0.0232 and 0.0184, respectively. This 
variation was about the same for all three injectors and seemed to be a result of manifold 
design. A satisfactory correlation was obtained for all stations when X was used as the 
characteristic dimension in Re, with C equal to 0. 0215. The correlation was 
applicable for pressures of 1.034 and 6.895 MN/m^ (150 to 1000 psia) and mixture ratios 
from 2. 57 to 8. 1 (28 to 11 percent hydrogen). The standard deviation of C was 0.0038. 

INTRODUCTION 

To effectively design convectively cooled nozzles, a detailed knowledge of the axial 
and circumferential variation of the hot-gas-side heat -transfer coefficients should be 



known. Reference 1 gave some of these details for a hydrogen-oxygen rocket. Since 
this work was published, numerous inquiries have been received for more detailed infor- 
mation on the axial and circumferential hot-gas -side heat-transfer coefficients, espe- 
cially in the straight section of the chamber and at the throat. Discoloration patterns on 
rocket chambers have also caused designers to worry about the effects that the injector 
and variations in mixture ratio have on the heat-transfer coefficient. 

To answer these questions, additional instrumentation was installed and further 
testing was conducted on two copper heat-sink thrust chambers. These thrust chambers 
had the same geometry as those of reference 1. Local heat-transfer rates were deter- 
mined at 20 locations (12 axial stations) for thrust chamber 1 and for 18 locations (five 
axial stations) for thrust chamber 2. Six circumferential instrumentation sites were in- 
corporated in the chamber and throat axial stations of thrust chamber 2. 

Three coaxial injectors were used. One had a copper faceplate and the other two 

had porous faceplates. One porous faceplate was designed for low chamber pressure 

2 

(nominal 2.068MN/m or 300 psia) and the other originally for high chamber pressure 

2 

(nominal 4.147 MN/m or 600 psia). The high-chamber -pressure porous faceplate in- 
jector was modified (its fuel elements had been damaged) and then used. This injector 
was run to determine the effect of misalined and damaged elements on the heat-transfer 
coefficients. 

o 

The thrust chambers were run at a chamber pressure of 2.068 MN/m (300 psia) 
over a range of mixture ratios from 2. 57 to 5. 67 (28 to 15 percent hydrogen). The data 
from reference 1 are also included in this report for correlation purposes. 

In the correlation procedure of reference 1, the station diameter was used as the 
characteristic dimension in the Reynolds number. In this report, the axial distance 
from the injector face is also used as the characteristic dimension. Comparisons based 
on these correlations are presented herein for the present rocket data, the rocket data of 
reference 1, and the heated -air data of reference 2. 


INSTRUMENTATION AND DATA RECORDING 

Copper rods, 0. 635 centimeter (0. 25 in.) in diameter and 6. 35 centimeters (2. 5 in.) 
long, were used to obtain the transient temperature data. The geometry, construction, 
installation, sealing, and pressurizing details of these rods are the same as those used 
in reference 1. Figure 1(a) shows the thermocouple stations on the copper rod. Four 
chromel-alumel thermocouples were used. Twenty rods at 12 axial stations were used 
on thrust chamber 1. Eighteen rods at five axial stations were used on thrust chamber 2. 
Six circumferential locations were incorporated in the chamber (station 2) and throat 
(station 8) axial stations for thrust chamber 2. A sketch of the chambers along with a 
table of the instrumentation coordinates is presented in figure 1(a). 
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(a) Instrumentation locations and sealing details. 
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(b) Copper heat-sink rocket thrust chamber with coaxial injector. ( All linear dimensions are in centimeters (in. ). ) 


Figure 1. - Instrumentation. 



Static pressure measurements were taken at all stations except at the 29. 49- and 
22. 21-centimeter (11.61- and 8. 744-in.) axial locations for thrust chamber 2. The mea- 
sured chamber pressure was used as the static pressure for these two stations, which 
are in the constant diameter position of the chamber. The theoretical static pressure 
was used at all stations for thrust chamber 2. 

A venturi and orifice were used to measure the oxygen and hydrogen flows, respec- 
tively. Temperatures and pressures were recorded on a digital recording system that 
had a sampling rate of 31 250 words per second with a block length of 125 words. Data 
parameters were fitted over a 25-word sampling interval with smoothing to eliminate 
60 -hertz noise and greatly diminish any random noise. A smooth curve was fit through 
25 readings of the data parameter, and then one reading, to be used in the terminal cal- 
culations, was made for all parameters at a common time. This reduced the amount of 
terminal calculations. Chamber pressure, which was sampled many times in each data 
block, was used as a triggering device for starting and stopping the calculations. These 
data were used to evaluate the performance as well as the heat transfer. The heat- 
transfer coefficients that were calculated from three thermocouples on each rod were 
averaged. The average was used in the correlation procedures. 


APPARATUS AND EXPERIMENTAL PROCEDURE 

In this investigation, two copper heat-sink, solid-wall thrust chambers were used to 
obtain transient temperature data. Figure 1(b) gives the geometry of the thrust cham- 
bers. 

Three coaxial injectors were used with liquid oxygen and gaseous hydrogen as the 
propellants. Each of the injectors had 234 injector elements uniformly spaced in a circu- 
lar pattern (0.4 element/cm or 2.6 elements/in. ). The injector faceplates were dish 
shaped. One faceplate (injector 1) was made of solid copper and the other two were made 
of porous Rigimesh stainless steel material. The Rigimesh was used to give the face- 
plate transpiration cooling (fig. 2). Injector 1 had oxygen tubes that had fins to center the 
oxygen element. No fins were used on the oxygen tubes that were used with Rigimesh 

faceplates. One porous faceplate injector (injector 2) was designed for low chamber pres- 

o 

sure (nominal 2.06 MN/m or 300 psia) and the other (injector 3) was designed for high 
chamber pressure (nominal 4.137 MN/m or 600 psia). Originally, the oxygen and hy- 
drogen hole sizes were the only difference between the two designs. The holes were 

2 2 
sized for a 689-kN/m (100-psi) pressure drop at chamber pressures of 2.068 MN/m 

(300 psia) and 4.137 MN/m (600 psia). Injector 3, which originally was the high cham- 
ber pressure design, was used in this heat -transfer study to find the effects of misaline- 
ment and concentricity of elements on the heat-transfer coefficients. The misalinement 


4 



I 



was caused by damage from many previous runs. The oxygen holes in this injector were 
bushed down to the same size as the low-chamber-pressure Rigimesh design so only the 
hydrogen feed holes were larger than injectors 1 and 2. Injector 3, which was a high- 
chamber -pressure design, was modified to run at low chamber pressures. Thus, the 
results from the injectors could be compared under similar operating conditions. 

The test runs were made on a test stand located at the Lewis Plum Brook Facility. 
Figure 3 shows the copper heat-sink thrust chamber installed in the test facility. Propel- 
lent valves for controlling gaseous hydrogen and liquid oxygen were positioned before the 
rim and opened to these fixed positions during the run to provide a step rise in chamber 
pressure. Once the nominal chamber pressure had been achieved, an automatic control 
took over and repositioned the valves to give exact values of chamber pressure and mix- 
ture ratios so that runs with different injectors could be repeated with identical flow con- 

9 

ditions. All runs were made at a chamber pressure of 2. 068 MN/m (300 psia) with 
mixture ratios varied from 2. 57 to 5. 67 (28 to 15 percent hydrogen). Full chamber pres- 
sure was achieved in 0.02 to 0.06 second. This step function in chamber pressure or 
driving temperatime made it possible to use a simpler mathematical model to obtain heat- 
transfer coefficients. 
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Figure 3. - Copper heat-sink rocket thrust chamber on test stand. 


CALCULATION PROCEDURE 

The calculation procedure used to find the heat -transfer coefficient was the constant 
h method as used in reference 1. (Symbols are given in the appendix.) 

The equation used for the solution of h from the one -dimensional semi-infinite slab 
as given in reference 3 is 



This equation was programmed so that an iteration process found the h that satis- 
fied the measured T for a given location x on the rod and time t. The initial condition 
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is T = T q at t = 0. The material properties (k,p, c) were evaluated at a temperature 
< T x=0 - T 0>/ 4 + T 0 < ref ' 4 >- 

The nondimensional heat-transfer parameters were computed by introducing the 
transport properties as a function of reference enthalpy H* (ref. 4) and static pressure 
P g where H* = H g + 0. 5 (H^ - H g ) +0.22 (Pr*)^^ (H c - H g ) . Numerous programs 
are presently in existence for the calculation of equilibrium compositions and other ther- 
modynamic properties of complex chemical systems. The programs of reference 5 for 
thermodynamic properties and of reference 6 for transport properties were modified and 
simplified for a gaseous -hydrogen - liquid-oxygen system by Frederic N. Goldberg of 
Lewis . 

Before the heat-transfer coefficient h can be computed, the driving temperature 
T aw must be determined throughout the nozzle; to determine T AW , the combustion 
temperature must be known. Combustion temperature is a function of combustion cham- 
ber pressure, percent fuel, and combustion efficiency. The program was further modi- 
fied to account for combustion efficiency and is described in reference 1. 

Two conditions were spelled out in reference 1: When the combustion efficiency was 
less than 1, an iteration was performed whereby combustion temperature was reduced 
until the measured weight flow and chamber pressure produced M = 1 at the throat; 
when the combustion efficiency was greater than 1, the theoretical combustion tempera- 
ture was used and the weight flow was adjusted to produce M = 1 at the throat. 

The average combustion efficiency obtained from all the test data was 98 percent of 
theoretical equilibrium. A 2-percent change in combustion efficiency results in a change 
of heat-transfer coefficient of about 2 percent. 

For thrust chamber 1, T A ^. was determined for each station in the convergence, 
throat, and divergence sections by using the measured static pressure ratio instead of 
area ratio. 


RESULTS AND DISCUSSION 
Axial Variation of Correlation Constant C 

A correlation of the form St* Pr*®‘^ = C was used in reference 1, where 

St* , Pr* , and Re^ are reference Stanton, Prandtl, and Reynolds numbers, respec- 
tively, and C is a constant for any one axial station. Figure 4(a) shows the variation in 
C with axial distance for 12 axial locations in thrust chamber 1 for various mixture 
ratios . 

The C variation with axial distance results for reference 1 and the mean values for 
the various mixture ratios for the present results are shown on figure 4(b). The values 
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(a) Correlation constant as function of axial distance for various mixture ratios. 
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(b) Correlation constant comparison of present data to that of reference 1. 



Figure 4. - Correlation constants. 


of C for reference 1 are the mean values for chamber pressures from 1.03 to 6.845 
MN/m^ (150 to 1000 psia) and for mixture ratios from 4.88 to 8.10 (17 to 11 percent hy- 
drogen). To show the spread of C, the 95-percent probability of the data from refer- 
ence 1 is also shown for these data. The data from the present tests are in good agree- 
ment with the results from reference 1 and show additional detail because more axial in- 
strumentation locations were used. The C values extrapolated to the injector face indi- 
cate higher values than would be inferred from reference 1 . 

At any one station, this equation does a good job of correlating the results over a 
range of fuel percentages using the mean value of C. In other words, the use of refer- 
ence enthalpy and static pressure for getting the transport properties in the correlating 
equation adequately accounts for variable properties. 


Circumferential Variation of the Correlation Constant C 

Table I shows the variation of h > C ’ T AW’ and Tq with circumferential location. 
Since the data were all taken by a transient technique, the reported values are not all for 
the same time for each injector; therefore, the absolute values of the wall temperature 
should not be used in comparing the results for the various injectors. The heat-transfer 
coefficient h or the correlating constants C can be compared. 

Figure 5(a) shows a plot of C as a function of circumferential location for the cham- 
ber (axial position, -22. 21 cm or -8.744 in.) and the injectors used with thrust cham- 
ber 2. The C was obtained from the correlation St* Re^'^ = C. Figure 5(b) 

shows a plot of the same variables for the throat position. Figure 5 shows that the data 
for all the injectors have about the same trends in amplitude and shape. The damaged 
element injector (modified high P porous, injector 3) gives the same variation in the 
constant C as the other injectors. The circumferential variations in h or C for the 
worst cases were 38 percent in the chamber and 27 percent at the throat. Regardless of 
how the thrust chamber was positioned with respect to the injector, a temperature sens- 
ing plug was always alined with an injector element. This was true even for the results of 
the modified high-chamber -pressure porous injector (3) where the nozzle was rotated 30° 
with respect to the injector. Therefore, the amplitude shown in figure 5(a) cannot be at- 
tributed to the misalinement of sensing plug and injector elements. 

Figure 5(a) shows an arrow indicating the gaseous hydrogen manifold inlet angular 
position. There is only one hydrogen inlet. This position coincides with the lowest C's 
for all injectors. Another region of low C's, 180° from this position, suggests that the 
manifolding might cause two hydrogen rich regions and two hydrogen starved regions in 
the chamber. 
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TABLE I. - PERFORMANCE AND HEAT TRANSFER 
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(a) Chamber correlation constants; axial distance from throat, -22. 21 centimeters {-8.744 in. ). 



Figure 5. -Correlation constants for various injectors. 
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Figure 6 shows a sketch of the injector showing the angular location of the hydrogen 
manifold inlet. The circumferential variation of heat input q, as determined from the in- 
strumented rods, is also shown on this figure for reading 24. Again, these variations 
are shown to be lowest at locations in line with and opposite to the inlet line of the hydro- 
gen manifold. 


Analysis of the Effect of Input Specifications on Calculations 
of the Heat-Transfer Coefficient 

The C's and q's presented in figures 5 and 6 were computed with driving tempera- 
tures that corresponded to the chamber pressures and percent fuels measured at the in- 
jector and listed in table I. Hereafter, this will be called input specification 1. 

To study this variation in more detail, two more input specifications will be studied 
using reading 24. These variations will only be studied at two angular positions corre- 
sponding to the regions of highest and lowest q's. 

Reading 24, for example, could be recomputed where the driving temperature T 
could be considered a function of the zone, that is, by allowing the percent fuel to vary. 
The highest wall temperature T^ (613.2 K or 1103.8° R) will give a heat-transfer coef- 
ficient h of 2.801 kW/(m 2 )(K) (0. 1371 Btu/(ft 2 ) (sec) (°R)). Assuming h constant cir- 
cumferentially , this h with the lowest wall temperature measured (510.7 K or 919.3° R) 
requires a driving temperature of 2541 K (4575° R) . This temperature corresponds to 
24 percent hydrogen. The h used here is the lowest h that can be used for both the 
highest and lowest wall temperatures measured. Any lower h would require a T^^, 
greater than stoichiometric to get the highest wall temperature. Hereafter, this type of 
calculation will be called input specification 2. 

Another way to handle the data would be to assume that a flat plate correlation of the 

n 7 -0 9 

form St* Pr* ‘ ‘ = 0.0295 Re* (hereafter called input specification 3) and that the 
correlating constant (0.0295) should be applicable at all circumferential positions at the 
station x = -22.21 centimeters (-8.744 in.). Then the measured values of q and T W 
could be used to determine the range of T^^- that would be necessary to satisfy these 
conditions. 

For reading 24 (table I) at an angular position of 285°, use of a T of 3340° K 
(6012° R) resulted in an h of 2.935 kW/(m 2 )(K) (0. 1437 Btu/(ft 2 )(sec)(°R)) and a corre- 
lation 

St* Pr* 0 ' 7 = 0.02838 Re* -0 ’ 2 
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Using a flat -plate correlation with axial length instead of diameter as the characteristic 
dimension results in 


St* p r * 0 - 7 = 0.02504 Re*" 0,2 

The heat -transfer coefficient h is relatively insensitive to changes in percent fuel. 
Therefore, the ratio of h's can be formed 

h _ 0.0295 
0.1437 0.02504 

to yield h = 0. 1693, which is consistent with a flat -plate correlation with a constant of 
0.0295. To satisfy a q of 8.008 MN/m 2 (705.6 Btu/(ft 2 )(sec)) and a wall temperature of 
613.2 K (1103.8° R), T aw would have to be 2928 K (5271° R) . With an h of 3.458 kW/ 

(m 2 ) (K) (0.1693 Btu/(ft 2 )(sec)(°R)), a T w of 510.7 K (919.3° R), and a q of 5.726 MW/ 
(m 2 )(sec) (504.5 Btu/(ft 2 )(sec)), would have to be 2166 K (3899° R). These temper- 

atures correspond to a mixture ratio of 4 and 2. 45 (20 and 29 percent hydrogen) . The 
mixture ratio supplied to the injector was 5.934 (15.64 percent hydrogen). 

The following table summarizes the three input specifications for calculating the cir- 
cumferential variations observed at station 2 in the chamber for reading 24. The calcu- 
lations illustrate that what is really measured in the rods is the heat flow rate q, and one 
has the problem of separating h and T^.^. Input specification 1 assumes a constant 
T AW at all circumferential locations. The is that computed from the measured 


Input specification 

Heat flow rate and 

Heat-transfer coefficient, h 

Adiabatic wall 

Percent 


wall temperature 



temperature, T A „ r 

fuel 


values 

kW 

Btu 


n vv 




(m 2 )(K) 

(ft 2 )(sec)(°R) 

K 

°R 


1. Constant T^ w , constant 

Highest 

2.935 

0.1437 

3340 

6012 

15.64 

percent fuel 

Lowest 

2.024 

.0991 

3340 

6012 

15.64 

2. Constant h circumferen- 

Highest 

2.801 

0.1371 

3340 

6012 

11.2 

tially, highest percent 

Lowest 

2.801 

.1371 

2540 

4575 

24 

fuel = stoichiometric 







3. Constant C circumfer- 

Highest 

3.458 

0.1693 

2928 

5271 

20 

entially, St* Pr*^‘^ - 

Lowest 

3.458 

.1693 

2166 

3899 

29 

0.0295 Re* -0 ' 2 

X 1 
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P c and percent fuel supplied to the injector. Specification 2 assumes a constant h at all 
circumferential locations and gives a variation of T^^ that goes along with a variation 
from stoichiometric (11.2 percent) to 24 percent hydrogen. Specification 3 assumes a 
constant C at all circumferential locations that goes along with a heat-transfer correla- 
tion of the form 


St* Pr* 0,7 = 0.0295 Re* -0,2 

X 


This results in a variation of T AW that goes along with a variation of mixture ratio of 
4 to 2. 45 (20 to 29 percent hydrogen. ) The measured mixture ratio at the injector was 
5.934 (15.64 percent hydrogen). 

The preceding calculations, although they do not pinpoint h exactly, do indicate 
rough limits that warn a designer that a good injection process and good manifolding must 
be used in order to avoid circumferential variations in heat transfer. Figure 5(b) shows 
that by the time the throat is reached, the pattern of C varying with angular position is 
no longer observed. These are the only two axial locations where detailed heat-transfer 
circumferential data were taken, and all figures showing axial heat-transfer detail indi- 
cate averages of the measurements at any one axial location. Since T^^ was not meas- 
ured, the driving temperatures used in the results for the rest of the report are those 
that go with the measured chamber pressure and percent fuel supplied to the injector 
(specification 1). 

To see how the circumferential variation would affect the wall temperatures in a liq- 
uid cooled thrust chamber design, a 15-percent variation in the hot-gas-side heat-tranfer 
coefficient was made. The thrust chamber used to analytically evaluate these wall tem- 
perature changes was a liquid -hydrogen cooled research chamber that had 347 stainless- 
steel coolant tubes with 25-millimeter (0. 010 in. ) wall thicknesses. The coolant tubes 

were designed to simulate a nuclear rocket; that is, all propellant went through the cool- 

2 

ant passages. The inlet conditions were 4.82-MN/m (700-psia) coolant pressure at a 
bulk temperature of 27.8 K (50° R). At a chamber pressure of 2.068 MN/m 2 (300 psia), 
this 15-percent variation caused the throat wall temperature to change by approximately 
111.1 K (200° R). The chamber and exit wall temperature changed approximately 
55. 56 K (100° R). These magnitudes of wall temperature variation could be important 
in a critical design. 


Correlation Using Reynolds Number Based on Axial Distance 

The previous correlations in which Re was based on the local diameter do a reason- 
able job of predicting h, provided that the local value of C is known (for a given station, 
the property evaluation at reference enthalpy handles various mixture ratios and chamber 
pressures). However, in a new engine geometry, the designer must determine the proper 
16 



axial distribution of the constant C. Van Glahn in reference 7, because of the complexity 
of nozzle heat transfer, gives several correlating equations in order to accommodate the 
several heat-transfer regimes as defined by Reynolds number and geometric locations. 

0 7 

In attempting to find a simpler correlation that could universally be used, St* Pr* 
data were plotted as a function of Reynolds number based on axial length from the injector 
face. In determining this correlation, the geometry of a high-contraction-ratio, large 
nozzle was assumed to more closely approximate a flat plate than a tube. Figure 7 shows 
plots of these data and also the data of reference 1 for the various axial stations. Fig- 
ure 7(c), which is the throat station, indicates that a slope other than -0.2 would better 
fit the data; however, the range of Reynolds number covered by the data is small. The 
middle 'v] symbols on this curve and that of figure 6 are the mean values of C for four 
runs at a chamber pressure of 2.068 MN/m (300 psia) and a mixture ratio of 5.67 
(15 percent hydrogen) for six circumferential stations. The upper and lower symbols 
are the average of the highest and lowest C's for these four runs. These limits are re- 
peated here to show the circumferential variation on a Re x basis. 

Figure 8 is a plot of the correlation constant C as a function of axial distance along 
the nozzle. The solid line connects C's computed from the equation 

C = St* Pr* 0-7 Re* 0,2 

X 

and the dashed line connects C T s computed from the equation 

C = St* Pr* 0 - 7 Reg 0 - 2 

When using x for the characteristic dimension in the Reynolds number, the variation in 
C is reduced by a factor of two. Figure 8 also indicates that the heat-transfer corre- 
lating constant C in the throat region is not much different from the average for the 
whole thrust chamber when put on a Re x basis. In other words, acceleration or other 
effects in the throat region for this thrust chamber do reduce the C used in predic- 
ting the heat-transfer coefficient, but only about 10 percent below the average C for the 
whole thrust chamber. Acceleration or other effects in the throat region of a nozzle when 
based on an Re using diameter as the characteristic dimension normally reduce C at the 

throat by approximately 40 percent below that value of C in the chamber (ref. 1) . 

0 7 

Figure 9(a) shows a plot of St* Pr* as a function of Re* for a chamber oressure 
2 x 

of 2. 117 MN/m (307 psia) and mixture ratio of 2.64 (27. 5 percent hydrogen) for all axial 

stations. Figure 9(b) is a plot of the same data as a function of Reynolds number based 

on the diameter. Figure 9 shows a typical correlating line for heating in tubes where 

0.026 is used as the constant. One sees that using x for the characteristic dimension 
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Figure 9 - Product of Stanton number and (Prandtl number) 0 - 7 as function of Reynolds number for various axial posi- 
tions. Reading 18; chamber pressure, 2. 119 MN/m 2 {307. 40 psia); percent hydrogen, 27.57; mixture ratio, 2.627. 
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spreads the data out along a line of -0.2 slope and that seemingly one line could be used 
for correlating all the stations. The slope would be far from -0.2 if one would attempt to 
correlate all stations with one line using Re^. 

Figure 10(a) is a plot of St* Pr*®' 7 as a function of Re x for all the data taken in 
reference 1 as well as the present data for all axial stations. A correlation of the form 

St* Pr* 0 ' 7 = C Re* -0 ' 2 

X 

correlates all the data for all stations for a C = 0.0215 with a la variation = 0.0038. 

The 95-percent probability lines of reference 1 are also shown on this plot. In figure 
10(b) the correlating equation and the 95-percent probability lines are replotted. The 'vl 
symbols on this plot are again the mean, the average highest C, and the average lowest 
C for four runs where data were taken at six circumferential locations. The spread is 
from 13. 6 to -9.8 percent at the throat and from 17.6 to -14. 1 percent in the chamber. 
The circumferential variation shown in this plot is for one running condition for the cham- 
ber and throat, respectively. The circumferential spread covers about half of the total 
spread of the correlation 


St* Pr* °- 7 = 0. 0215 Re* _0 ' 2 

X 

which covers all axial stations, all chamber pressures from 1.03 to approximately 6.895 
MN/m^ (150 to 1000 psia), and all mixture ratios from 2. 57 to 8.1 (28 to 11 percent hy- 
drogen). It also gives a designer some insight into how much of a safety margin would 
have to be allowed for in circumferential variations of heat transfer for a cooled thrust 
chamber design similar to this thrust chamber. 

In order to further check the reliability of this correlation, the heated air data of 

reference 2 were used. These data are especially interesting because the cooled chamber 

0 7 

length was varied. Figure 11 shows St* Pr* * as a functiDn of Re* for chamber 

lengths of 0, 15.24, 30.48, and 45.72 centimeters (0, 6, 12, and 18 in.), where x is the 

2 

distance from the chamber entrance. These runs were all heated air at 517. 1 kN/m 
(75 psia) and 883. 3 K (1500° R) . To increase the Reynolds number range, a run with a 
45.72-centimeter (18-in.) chamber length is shown where the total pressure was in- 
creased to 1.72 MN/m 2 (250 psia). The equation St* Pr* 0,7 = 0.0247 Re* - ® -2 seems 
to correlate the combined data for all stations with a la variation of 0.0047. Since the 
range of Reynolds number was small, the -0.2 slope customarily used for turbulent flow 
was retained. The 95-percent probability lines for these data are also shown on this plot. 
Figure 12 shows that using x instead of d for the Reynolds number makes the data fol- 
low a flat plate-like correlation. It also shows that acceleration effects do not reduce the 
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Figure 10. - Proouct of Stanton number and (Prandtl number) 0 * 7 as function of Reynolds number. 
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Figure 11. - Product of Stanton number and (Prandtl number)®- 7 as function of Reynolds number for 
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C to be used at the throat if an equation of the form St* Pr* ^ = C Re* ^ j S use d 
for predicting the heat-transfer coefficient. The authors of reference 2 used the turbu- 
lent boundary -layer analysis of reference 8 where 


p u 

• P I 


C p = K * 


C f* 


e e 


n 


K* = 



5Pr + 51n(Pr + 1) 



The factor K* is similar to the Prandtl number correction factor in the von Karman 
analogy. The coefficient Cjt is analogous to the wall friction coefficient Cj but with 
the momentum thickness replaced by the energy thickness. When n = 0, the Stanton 
number depends only on the thermal characteristic cp. Therefore, the preceding equation 
becomes 


h 

Pe u e 


P 



Using this equation, the authors of reference 2 predicted h for the 0-, 15. 24-, and 
45. 72 -centimeter (0-, 6-, and 18-in.) chamber lengths. Figure 13(a) is a plot of the 
ratio of h experimental to h predicted by the preceding equation for the various axial 
positions and the three different chamber lengths used in reference 2. A similar plot of 
the data of reference 2 is presented in figure 13(b) except that h experimental is divided 
by h predicted from a flat plate type of equation St* Pr*^‘ ^ = 0.0247 Re* ^ , where 
x was always measured from the beginning of the chamber. The same trends in the re- 
sults are apparent in figures 13(a) and (b) with the simpler flat -plate -like correlation 
doing a slightly better job. This flat -plate type of correlation needs to be tried on more 
nozzles of different geometries, propellants, etc. 
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Figure 13. - Ratio of experimental heat transfer coefficients to those predicted for various inlet lengths. 


CONCLUDING REMARKS 

Although this particular thrust chamber and injector design had circumferential var- 
iations in q, these effects were investigated to the extent that it was determined that a 
valid axial correlation could be developed. In the light of the findings of this report, a 
word of warming must be given to a designer. He must design for his nozzle geometry 
and his injector. He must allow for circumferential effects. Because the methods for 
generating these data by experiment are now proven methods and are not too complicated, 
it seems feasible that one of the steps in the development of any new rocket system should 
be a heat-transfer experiment rim early in the program perhaps in conjunction with injec- 
tor stability tests. This is true only for designs that push the local cooling capabilities of 
the system; otherwise, large safety margins can be used in design and these tests by- 
passed. Simple correlations like St* Pr*®'^ = C Re* ^ or St* Pr*^*^ = C Re^ 
can be used. The Re type of correlation is recommended for the chamber region near 
the injector, especially for large diameter chambers with high contraction ratios. 
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SUMMARY OF RESULTS 


An experimental investigation was conducted at NASA-Lewis to determine the axial 
and circumferential variations of heat -transfer coefficients in two rocket thrust cham- 
bers. Transient temperature measurements were made at 20 locations in one thrust 
chamber and at 18 locations in another thrust chamber. Six circumferential locations 
were used at the chamber and throat stations . The thrust chambers were operated over 
a range of mixture ratios from 2.57 to 5.67 (28 to 15 percent hydrogen) at a chamber 
pressure of 2.068 MN/m (300 psia). Three injectors were used. 

1. The data for all stations of this thrust chamber are correlated by the equation 
St* Pr*^'^ = 0.0215 Re* - ^'^ with a la variation of C equal to 0.0038, where St* , 
Pr* , and Re* are reference Stanton, Prandtl, and Reynolds numbers, respectively, 
and C is a constant. When Re^ is used in this equation, the variation of data was 
twice that using Re x , where d is the diameter and x is the axial distance from the in- 
jector . 

2. The maximum circumferential variation of C in the chamber and throat are 
38 and 27 percent, respectively. 

3. In the throat region C was equal to 0.019, this is about 10 percent below the av- 

0 7 

erage value of C (0.0215 for all stations) when based on the equation St* Pr* ‘ = 

C Re* -0 ’ 2 . 

4. The Jet Propulsion Laboratory heated air data for a nozzle gave a C = 0.0247 
with a la variation of 0.0047 for all stations where Re was used. 

X 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 22, 1971, 

122-29. 
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APPENDIX - SYMBOLS 


A 

A* 

c 

c f * 

c 

D 

H 

h aw 

h 

K* 


k 

L 

P 

Pr 


q 


Re 

Re 


d 

x 


S 

St 

T 

t aw 

t 


cross-sectional area 

throat cross-sectional area 

constant in equation St* Pr*®'^ = C Re* - ®'^ 

coefficient analogous to skin-friction coefficient, with momentum thickness de- 
pendence replaced by energy thickness nozzle diameter 

specific heat of material 

nozzle diameter 

enthalpy 

adiabatic wall enthalpy 
heat-transfer coefficient 
correction factor in the equation 


r 



L 


5Pr + ln(Pr + 1) - 14 + 



thermal coefficient of conductivity of material 

copper rod length 

pressure 

Prandtl number 

heat flow rate per unit area 

* d 

Reynolds number based on diameter, p*V//i* 
Reynolds number based on axial length, p*V x /p* 
entropy 

Stanton number 
temperature 

adiabatic wall temperature, P ) 

time 
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w total weight flow rate 

x distance normal to heat surface 

X axial distance from injector face 

X* axial distance from nozzle inlet 

9 momentum thickness 

p material density 

a standard deviation 

c p energy thickness 

Subscripts: 

c chamber or total 

d based on diameter 

e condition at free-stream edge of boundary layer 

min minimum 

w wall 

s static 

th throat 

x based on distance from injector 

0 zero burning time 

Superscript: 

* reference enthalpy condition 
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